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We study the problem of optimization over positive valued-operator measure to extract classical 
correlation in a bipartite quantum system. The proposed method is applied to binary states only. 
Moreover, to illustrate this method, an explicit example is studied in details. 

I. INTRODUCTION 

Quantum computing constitutes a rich research area of physics and computing at the same time. It is believed 
that, the expected power of a quantum computer is derived from genuine quantum resources. Entanglement, and 
correlations in general, are typical quantum resources. However, not all correlations have pure quantum nature. 
Generically, total correlations are "mixture" of classical and quantum correlations. An important issue is to know to 
what extent classical correlations are used in teleportation protocols and quantum algorithms. For example, if one is 
able to determine the classical part of correlations then by the optimal measurement he can extract some information 
under in classical form leaving the quantum state with less entropy. For this procedure to be useful it should be 
done while retaining the ability to regenerate the source state exactly from the classical measurement result and the 
post-measurement state of the quantum system. This has been studied in [1]. A possible application of this is to 
send the post-measurement state through a noiseless channel while sending the classical information through a more 
robust channel. 

Quantifying correlations implies measurement in most cases, which is a non trivial task for quantum states. A 
general measurement strategy is described by a positive operator-valued measure (POVM) that decomposes the unity 
in the Hilbert space of the particle under measurement. For the different measures of correlations (or information) one 
has to find the optimal POVM, which can be projective (orthogonal/nonorthogonal) or non-projective, and then the 
number of elements in the optimal POVM. POVM optimization is studied in different contexts in quantum information 

theory like the accessible information, the cost function. Fidelity, Davies' theorem [2] gives an upper bound for 

the maximum number of elements of an optimal POVM needed to attain the accessible information (or the mutual 
information) of an ensemble of states in a d-dimensional Hilbert space, this bound is d^ . The upper bound has been 
reduced to d{d+ l)/2 for states that are real [3]. There were attempts to reduce this bound to d elements, but Shor 
[4] has found an ensemble of states in a three dimensional Hilbert space that was optimized with six elements POVM. 
This eliminates the possibility of further reduction of the upper bound. Similar POVM optimization was made for 
the fidelity in [5], where for an ensemble of quantum states in 2-dimensional Hilbert the optimal POVM was found 
to have 3 elements, which fits the d{d -f l)/2 bound. In the present paper we present a POVM optimization for the 
classical correlations. The study leads to two important results: for arbitrary binary states projective POVM are 
found to be sufficient to optimize classical correlation which resembles the cases studied in [2,3], on the other hand 
the optimal POVM is found to have at most 4 (= 2^) elements with indications that the number can be reduced to 
3 (= 2(2 -I- 1)/2) for real states. Moreover, for specific states the optimal POVM is shown to be orthogonal. Another 
important property that is shown in the paper is the possibility of decomposing the mutual information contained in 
a given bipartite system into classical correlation and quantum discord. 

The paper is organized as follows: the definitions of POVM and classical correlation and their physical and math- 
ematical properties are given in the first two sections of the paper. Then we show the natural decomposition of the 
mutual information into classical correlation and quantum discord. The remaining part of the paper is devoted to 
the complete optimization procedure, including numerical simulation, of the classical correlation for a particular state 
illustrating the need for the quantum discord in order to match the mutual information. 
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II. POVM 



Let B = {Bi} be a POVM, then Bi should be a positive valued operator, to preserve the positivity of the mea- 
surement outcome probability. Each POVM "set" should decompose the unity , to get probability completeness, 
i.e. 

^B, = l. (1) 

For later use it should be mentioned that the set of all POVM's is convex, i.e. the segment joining two POVM's is a 
POVM. In simpler terms: if Bi and B2 are two POVM's and < A < 1 (a probability) then 

XBi + {l- X)B2 

is indeed a POVM. 



III. CLASSICAL CORRELATION 



Unfortunately, there is no unique measure of classical correlations. Different measures are found in the literature 
[6-11]. Two of the most relevant requirements that should be satisfied by any classical correlation are: (I) the 
classical correlations of product states {pab = PA ® Pb) should be zero. (II) Classical correlation should not be 
affected by local unitary transformations, which corresponds simply to a change of basis. Other important conditions 
are listed in [8,9]. Comparison between different measures is found in the literature. For example, property II has 
been investigated recently by B. Synak and M. Horodecki [12] for the measure proposed in [9], called the classical 
information deficit (A|^). It is found that does increase under local operations. Moreover it was shown [12] that 
is bounded above by the classical correlation measure proposed in [8]. The only case where A^ is monotone 
under local operation is when it coincides with the classical correlation of [8] . This comparison can be considered to 
be in favor of the Henderson and Vedral measure. However there is no evidence that the other measures of classical 
correlations violate some of the expected physical properties of a correlation. So one can choose a classical correlation 
measure out of the proposed measures. Our choice is the Henderson and Vedral measure. 

Given the bipartite state pab a possible measure of the classical correlation between subsystems A and B is ( [8]) 



Cb{pab) = max 



S{pa) - J2PiS{p\) 



(2) 



where pa = tvsipAB) is the reduced density matrix. The Von Neumann entropy is S{p) = — tr(plogp). S is a POVM 
and the sum over i runs over all its elements Bi. The conditional density matrix p\ is the density matrix of A after 
performing the measurement Bi on B: 

i _ trg {B.pab) , . 

tVAB{B,pAB) ■ ^ ' 

The probability of A being in the state p\ is pi = trAB{BiPAB)- 

The correlation measure Cb has a simple physical interpretation: if A and B are not correlated then the marginal 
entropy of A {S{pa)) and the residual entropy of A after a POVM measurement on B (J^Pi^iPli)) should coincide 
to give Cb = 0, since for uncorrelated system AB, A is not affected by a POVM measurement on B. Moreover, note 
that 



hence for a given POVM B the combination 



PA = ^PiPA 



S{pA)-J2piS{p\) 



is closely related to the entropy defect defined by Levitin (see [13] and references therein). So the classical correlation 
C B can be seen as the maximum average decrease in the entropy of the system A when a state (after a measurement 
Bi is preformed on B) is specified compared with the situation when only the mixture of states pA is known. 
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The correlation Ca{pab) is obtained from Cs by making the replacement A ^ B in the above formulas. It is 
evident that the measure Ca,b are not manifestly symmetric under the exchange of the roles of A and B. Whereas, 
one would expect that the classical correlation is a measure of how strongly the two subsystems are correlated no 
matter which subsystem is used to extract such a correlation. Hence, until a formal proof of an explicit symmetry of 
this measure, if it exists, the above proposed measure can not be considered as a universal measure of the existing 
classical correlations. 



IV. POVM CHOICE AND CONVEXITY 

In what follows we consider binary states exclusively. Extension to higher dimension is not straight forward. 



A. Convexity 

An important property that helps in reducing the range of exploration in the set of all POVM's is the concavity 
[14] of the classical correlation measure. If one is able to show that the classical correlation is convex then (see 
[2]). for binary states, one can consider POVM's of rank one elements that can be taken to be proportional to the 
one-dimensional projectors. 

Proof: 

For binary states the Von Neumann entropy is concave i.e. 

5(Api + (1 - X)p2) > XS{pi) + (1 - X)S{p2) (4) 

for all < A < 1. 

To study the concavity of the classical measure it suffices to study the variation of the functional F{pA,B) defined 

as 

F{pa,B) = -pbS{p^), (5) 

where the compact notation: ps = tiAB {BpAs) and = ^r^ {BpAs) /pb is used for all elements Bi of the POVM 
B. 

Let C and V be two POVM's and let Q be their combination 

g = XC+il-X)V, (6) 

for < A < 1. is a POVM since, as mentioned before, the set of all POVM's is convex. 
It is straight forward to show that 

= ApC + (1 _ a)p2 (7) 

pi = A4pA + (l-A)4/'A. (8) 
Pa Pa 

Using these equations and the fact that the entropy is a concave function, it can be shown that 

FipA.G) < XF{pA,C) + (1 - X)FipA,V). (9) 

Therefore, _F is a convex over the convex set of all POVM's. Hence the maximum of F occurs for an extremal POVM 
[15]. For binary states it was explicitly shown recently [16] that extreme POVM's arc of rank one^. Hence, to optimize 
F and consequently the classical correlation, since S{pa) is independent of the POVM performed, one can consider 
the special class of POVM with rank one elements. Hence we are left with projective measurements only^. 
An important consequence of this result is presented in the next section. 



'^For higher dimensional Hilbert spaces there are higher rank extreme POVM's [16]. 

^It should be mentioned that for generic cases one can use Naimark theorem to get projective POVM only without being 
restrictive, but this is done at the cost of extending the original Hilbert space. 



3 



B. Quantum discord and mutual information 



As mentioned in [12] the quantum mutual information I{A : S) of a bipartite system AB can be decomposed into 
an information deficit (or work deficit) A and the classical information deficit Ad^: 

7 = Aei + A . (10) 

So, naturally one should be concerned with an analogous decomposition of the mutual information when using a 
different measure of classical correlation. 

By analyzing different quantum states, the authors of [8] found that the estimated classical correlation and the 
relative entropy of entanglement {Ere) do not add up to give the Von Neumann mutual information between the two 
subsystems, i.e. 

(CB(pAB))opti^i,ed + Ere < I{PA:b) • (11) 

It was argued that either the mutual information is not the best quantity to measure the total correlations, or that 

a more elaborated choice of the POVM may saturate the total correlations. By considering all possible POVM for a 
given state, we show that the most optimal POVM can not saturate the total correlations. This implies that: either 
the mutual information is not a good measure of total correlations, which is highly improbable, or one has to change 
the measure Cb of classical correlation. The third alternative, which renders the definition of classical correlations 
compatible with the mutual information, is the possibility of having a different definition of quantum correlations 
which is different from Ejie- A possible candidate is the quantum discord, which was first defined in [17]. It is 
the result of the difference between classical and quantum conditional entropies. In contrast to classical conditional 
entropy, quantum conditional entropy is a measurement dependent quantity (see further [18]). The quantum discord 
is defined as [17] 

S{A : B) = I{A : B) - J{A : -B){nf } • (12) 

J is the information gained about B as a, result of the set of measurements {nf }: 

J{A : S){nf } = S{A) - S{A\{nf}) = S{A) - ^ *Pi5(V^) . (13) 

where *pi and *p\ were defined originally as the special case of pi and p\ when the choice of the set of measurement is 
restricted to one dimensional projectors 11^. The above definition of the quantum discord reflects clearly the inherited 
measurement dependence of the quantum conditional entropy. For the quantum discord to be zero one has to find at 
least one measurement for which it is zero. Therefore the minimum of the quantum discord is the relevant quantum 
correlation. By focusing on perfect measurements of B, defined by a set of one dimensional projectors, one can easily 
check that the set of measurements that minimizes the quantum discord, i.e. maximizes J, is exactly the same POVM 
set that optimizes the classical correlations for binary states. This follows from the definition of both quantities, and 
the result of the previous section on the optimization of classical correlation using projective measurements only. 
Hence 



Max{J{A : B)^uf}) = Cb{pab) ■ (14) 

Therefore 

I(A ■.B) = Cb+ min 5(A : B) , (15) 
{nf} 

i.e. for binary states the classical correlation and the quantum discord add up to give the mutual information. 



^The previously mentioned classical information deficit is a one-way quantity i.e. for a system with two parties A and B 

we allow communication from Ato B only or vice versa. This asymmetric definition is analogous to the seemingly asymmetric 
property of the classical correlations where one has to perform measurements on one of the two parties. 
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C. What is next? 



Limiting the domain of exploration in the set of all POVM's to projective POVM leads to important simplification. 
But the game is not over. What would be the optimal projective POVM? How many elements there arc in the optimal 
POVM? Corollary 1 of [16] implies that POVM's with 5 or more elements acting on a two dimensional Hilbert space 
are not extreme, and hence cannot optimize the classical correlation*. So, we are left with POVM's having 2, 3 or 
four elements. In the next section we illustrate, through an example, the steps that can be followed to complete the 
optimization procedure. 

V. AN EXAMPLE 

Consider a system AB in a state pab with probability p to be in a state pi and (1 — p) to be in another state p2- 
A possible example is the following state studied in [8] 

p\0){0\ ® |0)(0| + {l-p)\+){+\ ® |+)(+| , (16) 

where |+) = (|0) + |l))/\/2. To enlarge the "library" of computed classical correlations we will not use the above 
state but rather a slightly modified version of it, namely 

PAB = P|l)(l| ® |1)(1| + (1 - ® |+)(+| . (17) 

The classical correlation of the state in eq. (16) was calculated in [8] using optimization over all orthogonal measure- 
ments, which was not a justified restriction. 

For the state given in eq. (17), it is useful to use the Pauli matrices decomposition. Hence the measurement 
operator Bi can be written as: 

B. = (i±f^) . (18) 

The projection direction is set by the unit vector bi. We should further require the completeness relation (eq. (1)): 

Y,h = 0, ^|6,| = 2. (19) 

The state pab in eq. (17) becomes 

PAB = ||1)(1| ® (1 - a,) + ^|+)(+| ® (1 +a,) . (20) 
It is straight forward to show that the density p\ of subsystem A, after the measurement on B, is 

p^ = a,|l)(l| + (l-aOI+)(+l (21) 

where 

P(l - k ■ e^) 

cti = . (22) 

p{l-bi-e^) + {l-p){l + bi-e^) 
Hence the entropy of the subsystem A after the measurement Bi of the subsystem B is 

5(p^)=Mc^f^) (23) 
where as usual h{x) = — a;log(a;) — (1 — x) log(l — x) is the Shannon's entropy function and 



*This is analogous to the bound found in [2]. 
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4^'^'^ = 1(1 ± ^af + (l-a02) . (24) 

are the eigenvalues of p\. 

The probabiUty Pi to have the subsystem A in the state p\ is 

Pi = tvAB {BiPAB) = ^ (p(l - k ■ e,) + {l-p){l + bi- e^)) . (25) 

Note that, the same procedure can be appUcd for the state given in eq. (16). One has to replace Gz by —Gz in 
eq. (20). This change of sign will propagate to all the subsequent equations. 

A. Choice of basis 

The previous section seems to underline a choice of the basis for projection. One can argue that the component 
{bi)y does not show up in the conditional density matrix p^, hence we can set it to zero in the projection operator. 
In other words measurement can be performed in the xz plane. In general, this is not correct it was shown in [16] 
(see section III) that an extreme POVM with four outcomes (n = 4) cannot have its four vectors hi coplanar. On the 
other hand, it is easy to show that the entropy 

is maximum, i.e. it is equal to the entropy of the reduced density matrix pA 

Y,PiS{p\) = S{pA) 

for a projective measurement in the y direction, since the state of B is completely random in that direction. However 
for classical correlations the minimum of X^Pj-^CPa) i'' needed not its maximum. Therefore, it is justified to use 
measurements such that {hi)y = 0. Hence {bi)z = cos 9i and {bi)x = sin^,. But, with this assumption the n = 4 
extreme POVM is excluded. This should be investigated further to see whether it is a general property of real states, 
to which the present state belongs, i.e. whether optimal POVM's for classical correlation of real binary states have at 
most 3 outcomes in analogy with the <i(d+l)/2 bound found in [3]. In fact, the proof of [3] (see Lemma 5 and theorem 
1) relies on the convexity of the accessible information and some of it scaling properties under group transformation 
this seems to be equally valid for the classical correlation. This will be studied in a future work. But until such a 
rigorous proof is given the present " exclusion" of the n = 4 case should be considered as a reasonable assumption. 

B. Optimization 

1. Method 1: Lagrange multipliers 

The determination of the classical correlation Cb requires finding a POVM that minimizes X^Pj'S'(/9^), given the 
constraint equations (19). Hence a priori we have to explore all possible POVM's and then picking up the one that 
minimizes the considered sum. The result of section IV A has limited the domain of exploration to projective POVM 
only. The optimization problem can be done using the Lagrange multipliers method. It is more convenient to recast 
the problem as 

min ^Pi5(p^) =^rif{ei) 
subject to rj — 2 = 
Vi cos 6i = 

^riSin6'i = 0, (26) 

where = and 

= ^ b(l - cos9i) + {l-p){l + smei)]hiu;l'\9i)) . 
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Let Ai 2.3 be the Lagrange multipliers associated with the above constraints. Hence the minimization problem is 
equivalent to find the values of ^i, r,, Ai^2,3 satisfying 

f{0^) + Xi+X2Cos0i + X3sm0i = O (27) 

ri(^^^-A2sin^i + A3Cos^i) = (28) 

^r, = 2 (29) 

^ricos6'i = (30) 

^rismei = 0, (31) 

This optimization should be done for n = 2, 3. 



2. Orthogonal projective measurements 



For two-elements POVM (n = 2) wc recover the orthogonal projective measurement''^, since 61+62 = 0. Hence 
O2 = 01 + The best measurement, i.e. the optimal measurement among the n = 2 POVM's only, is obtained 
when f{6i) = /(^2) = f{0i + ""). Therefore, the best set of measurement is a two-shot measurements with opposite 
directions. This should be solved numerically for each p. The result is plotted in figure 1. Our result matches well with 
the result of [8] . This implies that although our state is slightly different from the state they have studied but the two 
states have the same classical correlations. However, as mentioned before, this is just the orthogonal measurement, 
non orthogonal measurements should be considered which was not done in [8]. 



3. Method 2: Monte Carlo simulation n = 3 



Finding the minimum by the method presented above (method 1) for n = 3 seems to be a non trivial task, however, 
it is possible to use a different method like the Monte Carlo Simulations (MCS) method. We generate random 
events that correspond to the set of the 3 angles {0}^^^ while from eq. 26 the three ri's are no longer independent 
variables. Thus our strategy is to compute the classical correlations for each randomly generated event. To each event 
a numerical value is obtained for the classical correlation, for a given p, this is is represented by the large number of 
points (forming long continuous strips) in Fig. 1. In Fig. 1 we plot the simulation outputs for 10^ events. Clearly 
the classical correlation predicted for the best orthogonal projective measurement case constitutes the upper limit 
of all the results obtained for the random events. Hence the optimal POVM, among all possible POVM's, is the 
two-outcomes orthogonal projective POVM. This result can be linked to the remark found in the conclusion of [13] 
that one should dare conjecture the following: orthogonal measiircmcnt is optimal whenever the number of states 
(the two states |1)(1| ® |1)(1| and |+)(+| ® |+)(+| in the present case) does not exceed the dimensionality of the state 
space. 



4. Mutual information 



To understand the role of the classical correlations and its compatibility with the mutual information, the mutual 
information is evaluated using the formula 

I{A : B) = S{pa) + S{pb) - S{pab) ■ (32) 



'""If one takes 61 = 62 = 1, it is straight forward to chock that the optimal POVM measurement is indeed projective, i.e. 
Bi = \^i){tpi\ and B2 = |V'2)('!/'2l, with the bases of measurement parameterized by 9 as 

{cos(|)|0) + sin(|)|l),sin(|)|0)-cos(^)|l). 
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For the present state of the system AB the mutual information can be easily calculated to obtain: 



I{A : B) = 2h{ ) - h{ ) . (33) 

I{A : B) is plotted in Fig. 1 besides the classical correlation obtained by the optimal orthogonal projective mea- 
surement and the one generated by MCS. Since pab is separable, then its relative entropy of entanglement is zero. 
Therefore, from our general POVM analysis, and as shown in Fig. 1, no POVM can lead to classical correlation that 
saturates, when added to the relative entropy of entanglement, the mutual information between the two subsystems. 
This confirms the necessity to use the quantum discord as a quantum counter part for the classical correlation used 
in this paper rather than the relative entropy of entanglement. 




FIG. 1. Wc plot as a function of p: the classical correlation C'b (lower curve) evaluated using the best orthogonal POVM, 
and the mutual information I {A : B) (higher curve). The vertical strips are sets of large numbers of dots representing the 
obtained value of the classical correlation for the different randomly generated events. 



VI. CONCLUSIONS 



By considering a generic POVM we show that classical correlations of binary states are optimized via projective 
POVM. It is found that the classical correlation and the quantum discord add up to give the mutual quantum 
information. 

This work should be considered as a first step towards a generalization of the POVM optimization for states in 
higher dimensions and to see whether rank one POVM elements are still the optimal choice for classical correlation. 
Moreover, there is at present strong evidence on the possibility of an experimental realization of non projective POVM 
using optical devices as was recently proposed in [19]. This gives a new dimension to POVM optimization. 

This work was performed as part of the research program of the Stichting voor Fundamcntccl Onderzoek der Materie 
(FOM) with financial support from the Nederlandse Organisatie voor Wetenschappelijk Onderzoek . 
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